Non-equilibrium Plasmons in a Quantum Wire Single Electron Transistor 
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We analyze a single electron transistor composed of two semi-infinite one dimensional quantum 
wires and a relatively short segment between them. We describe each wire section by a Luttinger 
model, and treat tunneling events in the sequential approximation when the system's dynamics can 
be described by a master equation. We show that the steady state occupation probabilities in the 
strongly interacting regime depend only on the energies of the states and follow a universal form 
that depends on the source-drain voltage and the interaction strength. 

PACS numbers: 71.10.Pm,73.23.Hk,73.63.-b 
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Strongly interacting one dimensional (ID) electron sys- 
tems that can be described by Tomonaga-Luttinger (TL) 
model 13] are under great interest due to recent experi- 
mental observations especially in the systems of single- 
walled carbon nanotubes (SWNT) HHllli. The Lut- 
tinger Liquid (LL), which is generahzation of the TL 
model, is characterized by the absence of fermionic quasi- 
particles; instead, the elementary excitations take the 
form of bosonic charge and spin density waves . 

Motivated by recent experiment by Postma et al. 0, 
we theoretically investigate a ID single electron transis- 
tor (SET). Schematic description of the system is that a 
finite wire segment, later called a quantum dot, is weakly 
coupled to two long wires as depicted in Fig. ^ Two 
point-like impurities or other defects in a long nanotube 
may result in equivalent physics. We model the system 
as a finite LL segment and two semi-infinite LL leads. 
The leads are connected to electron reservoirs which keep 
them in internal equihbria. The chemical potentials of 
the leads are controlled by source-drain voltage {Vsd), 
and we assume that tunneling barriers between leads and 
dot are high enough that sequential tunneling is the dom- 
inant charge transport mechanism. Closely related stud- 
ies based on a different set of approximations were re- 
cently published by Braggio and co-workers 8, 9]. If the 
tunneling barriers are lower, or if sequential tunneling is 
blocked (Coulomb blockade), correlations between tun- 
neling events across the left and right junctions become 
important. Transport in this regime has been discussed 
by Thorwart et al. jlol |. 
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FIG. 1: Model system. Two long wires are adiabatically con- 
nected to reservoirs and a short wire is connected to the two 
by weak coupling. 



In this Letter, we analyze the steady-state occupation 
probabilities of the many-body states on the dot. We 
will show that in the strongly interacting regime they are 
characterized by a universal distribution that depends on 
the applied voltage and the interaction strength. 

Carbon nanotubes have four transport channels in- 
cluding spin degeneracy 0, Q • The electron transport 
through quantum wires (QWs) is determined mostly by 
charge effects so it is reasonable to consider, as a first 
approximation the contribution of the total charge 
sector to SET characteristics, and postpone the discus- 
sion of the role of the other channels to a later date. The 
effective Hamiltonian H = Hq + Ht f^or this branch is 
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Here Hq is the bosonized LL Hamiltonian and the tun- 
neling Hamiltonian Ht describes bare electrons tunneling 
through barriers at the points Xl and Xfj, where indices 
i=l,2 and (t=± identify bare electrons. The parameter 
g indicates interaction strength, g=l for non-interacting 
electrons and g—>0 for strong repulsive interaction. The 
term Ec=^Tjf^{N — Nq)^ is on-dot zero-mode contribu- 
tion including the Coulomb charging energy, where Nq 
is a dimensionless gate voltage control parameter in the 
range 7VGe[0,l]. Leads are assumed to be semi-infinite 
with continuous spectra; the dot has finite size Lu and 
discrete energy spectrum with the plasmon level spacing 
AE — ttHvp / gL]j. 

We describe the system's dynamics using a master 
equation within a sequential tunneling regime where un- 
correlated single electron hoppings are dominant trans- 
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port mechanism [lllIT^. 
{d/dt)P{N,{n}) 

N'=N±1 {«'} 

-P{N,{n})T{N,{n}^ N',{n'})], 



(2) 



where the variables {N, {n}) represent a state of the dot 
with iV excess electrons and a set of the occupation num- 
bers {n} = (ni,n2, ■ ■ ■ ,n.„j, • • • ) of the plasmon states 
with momenta pm — m ■ ttH/Ld, P are the occupation 
probabilities of the dot-states [N, {n}), and T — -t-Ffl; 
are the transition rates due to tunneling through junc- 
tions at Xl and Xr. The lead degrees of freedom were 
integrated out since the leads are locally in equilibria. 

If the transition rates are known, it is straightfor- 
ward to determine steady-state occupation probabilities 
by solving master equation (|2J) for {d/dt)P{N, {n}) — 
with the constraint J2n {«} ^(^' {"■}) = 1- 

The transition rates r^/fl are evaluated within the 
golden rule approximation with respect to the tunneling 
Hamiltonian Ht using the standard bosonization tech- 
nique 0. As a function of the states on the dot and 
their energies the rates are given by 

TL/R{N,{n} N\{n'}) = 87r|ti/«|V?i 
xj{Ed{N', {n'}) - Ed{N, {n}) T {N' - N)eVL/R) 
x\{N', {n'}\i;\y6N',N+i + ^dSn',n-i\N, {n})\^ 

(3) 

where "— sign corresponds to index "L/R" and the 
voltage drops across junctions are Vl/r = j^^pl^^Vgd 
for the tunneling resistances Rl/r oc |iL/7?l^^ when dc 
source drain voltage Vsd is supplied. Lead contributions 
are given by a spectral function 7(e) |^ ^ 
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Here a = {g ^ — l)/4, appropriate for an end-contacted 
four-channel wire such as SWNT /3 = l/fcsT is the 
inverse temperature in the leads, and Ag = Ag^/^^~B^ 
where A is a short wavelength cut-off. At zero tem- 
perature, 7(e) is proportional to a power of energy, 
lim7-^o7(e) 6(— e)|e|". The dot energy Eo is ob- 
tained by solving the eigenvalue problem of the Hamilto- 
nian and is given by a sum of the zero mode energy 
EciN) = {T:hvF/2g^LD) ■ {N - Nq)^ and the plasmonic 
excitation energy E{{n}) — {nhvp / gLo) ■ J2m ' ^"i' 
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The zero-mode overlap is unity for A^' = A'^ ± 1 and van- 
ishes otherwise, and the overlaps of the plasmonic states 



can be written in terms of Laguerre polynomials as 

l(K}|^LlM>P = ^ff^' 
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where = min(nj„,nm) and n™'' = max(n^,n„i). 

Notice that at low energies only the first few occupations 
n„j and n'^ in the product in Eq. ^ differ from zero, 
yielding non-unity factors to the transition rate. 

With the analytically calculated transition rates (Eqs. 
we solve the master equation numerically and 
determine occupation probabilities of the many-body 
states. Hereafter we will present numerical results for the 
case with symmetric barriers, the dot size oi = 50nm 
and the Fermi velocity of vp = 8 x lO^m/s 

Fig. [3 shows (A) the occupation probabilities P{{n}) 
vs. the energies E{{n}) of the many-body states of 
(iV=0,{n}) and (B) the averaged probabihties = 
{P{{n}))n of degenerate states with energies E{n). The 
degeneracy D{n) of E{n) = n ■ ttHvf /gLo is the number 
of the many-body occupation configurations {n} that sat- 
isfy E{{n}) — E{n), i.e., J2m ' ^™ ~ 'which asymp- 
totically follows the Hardy-Ramanujan formula 
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As seen in Fig. |2IA), there is a dramatic change of 
the occupation probability distribution between a non- 
interacting {g—1) and a strongly interacting {g < 0.24) 
quantum dot. In the non-interacting case the occupa- 
tion probabilities of degenerate many-body states vary 
widely as is natural since electrons tunnel independently 
of each other (apart from the Pauli principle); typically 
states with many long wavelength plasmons are more 
likely than states with few short wavelength excitations. 

In contrast, in the strongly interaction regime all de- 
generate states are almost equally likely (variation of 
about 10% only, see the inset of Fig. IHB)) and the oc- 
cupation probability decreases nearly exponentially with 
energy. The suppression of variations originates from in- 
teractions in the central island [l^ and can be under- 
stood as a consequence of the separation of zero-modes 
and bosonic excitations in the strongly interacting limit 
even for a small system. The inclusion of spin effects 
in the central segment (spin is already included in the 
leads) introduces non-interacting degrees of freedom that 
results in increased fluctuations which can be suppressed 
by a magnetic field. In Fig. [2lB), we show that the 
average occupation probabilities P{n) = {P{{n}))n de- 
cay nearly exponentially even in the weakly interacting 
regime, g > 0.4, where the probabilities of degenerate 
states fluctuate strongly. 

The nearly exponential decay of the average probabil- 
ities suggests a Gibbs-like distribution. However, there 
is no obvious reason why this non-equilibrium distribu- 
tion should follow an equilibrium form. We obtain an 



3 



10" 



10" 



■10 



10" 





(A) 


g=0.24 




e 

8 


• 




1 


1 




1 


10^ 


g=i 


' « , 


0.2 0.4 p 



0.5 1 1.5 

E({n}) [eV] 

.8 g=,l-0O, 




FIG. 2: (A) The occupation probabilities P{{n}) = 
P{N=0, {n}) as a function of the energies E{N=0, {n}) of 
the many-body states on the dot for g = 0.24 (SWNT) and 
for g = 1 (non- interacting) in the inset. The degeneracy 
of the energy state E{n) is D{n) ~ exp{^v^/2n/3)/{4:\/^n), 
where n = X^m ' dimensionless energy (cir- 

cles (o) represent individual many-body states). (B) Aver- 
aged probabilities -P(n) — {P{{n}))n as a function of the 
energy levels E{n). Solid lines are analytic approximations 

P{n) oc n ^ "sd where n^d — eVsd ■ {irhvp / qLd) ■ Inset 
shows the relative standard deviations of the probabilities, 
{(5P({n})),i/P(n), as a function of the dimensionless energy 
n — E{n) ■ {-Khvp / gLD)~^ for g-values ranging from g — 1 
to g — 0.12. Implicit parameters are T = 50mK,Vsd ~ 
g-^ X 0.22V, and Ng=0.5. 



analytic approximation to P{n) at T=Q by setting the 
on-dot transition elements in © to unity and consid- 
ering the scattering-in and scattering-out processes for 
a particular many-body state {n}. In this approxima- 
tion, all degenerate states have equal probabilities and 
the total scattering rates at a specific energy are given 
by the product of the lead contribution which fol- 
lows a power- law, and the degeneracy Q). The master 
equation now reads 

d 

— P(n) = '^[P{n')D{n')T{n' n)-P{n)D{n')r(n r, 

n' 

(8) 

where the total charge labels N have been suppressed. 



By solving the master equation using saddle point inte- 
gral approximations, we obtain a differential equation for 
ln[P(n)] and find that for a large E{n), P{n) is 
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where Ugd = eVsd • {irhvp / gLjj)^^ is the dimensionless 
bias voltage and Z is the partition function. This distri- 
bution is universal in the sense that it does not depend 
on the precise form of the Hamiltonian of the dot as long 
as the excitation spectrum is bosonic and has linear dis- 
persion. Over limited energy ranges, it resembles a Gibbs 
distribution with the effective temperature 
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In the range of our interest 
SWNT device 0) and ti 



3 a + 1 1 + log n 
(10) 

g < 0.24 {g ~ 0.24 for a 
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10, this is approximately 
C{n) ■ g ■ eV,d/2 (11) 



where 1 < C{n) < 2. 

As an application of the newly found distribution func- 
tion, we investigate the current and the differential con- 
ductance G=dI/dV of the SWNT-SET with g = 0.24, as 
a function of gate voltage. Current / is 

^ = [^(^' {n}^N + 1, {n'}) 

N,{n},{n'} 

- P{N + 1, {n'})rL{N + 1, {n'}) ^ iV, {n}) . 

(12) 

Fig. O shows the /- characteristics: (A) off-Coulomb 
blockade with Nc=0.5 and (B) a blockade case with 
A^G=0.35 (/ is in units proportional to I^L/flP)- 

At zero temperature, there are conductance peaks 
whenever voltage drop across a junction equals differ- 
ence between two quantized energy levels. At a finite 
temperature only peaks that overcome thermal broaden- 
ing survive. For A^g'=0.5 (Fig. OlA)), electron numbers 
iV=0 and A^=l have degenerate ground states which re- 
sults in lifting of Coulomb blockade at zero bias voltage 
and yields equidistant conductance peaks separated by 
2AE=2{TThvF/ gLo)- In Fig. OjB), there is no trans- 
port at low bias voltage due to Coulomb blockade. The 
main conductance peaks indicate transitions involving 
the ground state. The side peaks (a) in Fig. OJA) and 
(a) to (c) in Fig. are for transitions between excited 
states 0, ^3 • The satellite peaks are smeared out at a 
high temperature as shown in the T=77K plot because 
they result from transitions between excited states. 

The inset in Fig. |3IA) shows conductance at zero 
bias voltage is proportional to a power of temperature, 
Gv^Q+iT) cx T"*, a* = a - 1 when ksTKAE and it is 
eventually expected to reach 2a when ksT^AE where 
the dot displays continuous spectrum. It can be general- 



ized to G 



(T) oc r"-i when ksTKAE. 



The G- V analysis agrees with Braggio et al. [a ISlI > 
noticing that our system has four channel exponent. 



4 




„xio-" 



^ G,50mK 

G,77K 

l,50mK 
I,77K 


,'(B) 

/ 






(c) / 














h 







°0 0.2 0.4 0.6 0.8 1 

vjv] 



FIG. 3; Current 7 and differential conductance dl /dV as a 
function of source-drain voltage Vsd for g = 0.24. (A) Off- 
Coulomb blockade with Ng=G.^ and (B) a Coulomb block- 
ade with 7Vg=0.35, for which the conductance vanishes up 
to eVad = 0.3 X Trhvp/g^ Ld at zero temperature. I is in units 
proportional to I^l/hP. Inset in (A) shows power-law of con- 
ductance on temperature at zero bias voltage where solid line 
is G{T) = G{To)iT/To)°'~^ for ksTKAE. Implicit parame- 
ters are same as in Fig. |2| 



a={g~^ — l)/A and spin effects in the dot are not con- 
sidered Since the results by Braggio et al. were 
obtained in a different approximation for the occupation 
probabihties, this comparison also shows that a simple 
/- V measurement is rather insensitive to the on-dot dis- 
tribution function. However, the distribution function 
is expected to have a more pronounced impact on sen- 
sitive probes such as noise measurements, which will be 
considered at elsewhere. 



In conclusion, we have analyzed a quantum wire single 
electron transistor where each wire segment is described 
by a Tomonaga-Luttinger model. We have shown that in 
the sequential tunnehng regime the steady-state occupa- 
tion probabilities of the many-body states at the central 
dot follow a universal distribution that is determined by 
a competition between energy dependent tunneling rates 
and the degeneracy of the on-dot many-body quantum 
states. We have applied the non-equilibrium distribu- 
tion to calculate the current-voltage characteristics of a 
SWNT-SET, and shown that the I-V measurement is 
rather insensitive to the distribution function. 
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